Spontaneous pattern formation in driven nonlinear lattices 
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We demonstrate the spontaneous formation of spatial patterns in a damped, ac-driven cubic Klein- 
Gordon lattice. These patterns are composed of arrays of intrinsic localized modes characteristic 
for nonlinear lattices. We analyze the modulation instability leading to this spontaneous pattern 
formation. Our calculation of the modulational instability is applicable in one and two-dimensional 
lattices, however in the analyses of the emerging patterns we concentrate particularly on the two- 
dimensional case. 



Complex spatial patterns are often observed in systems 
driven away from equilibrium !§. Typically, the pat- 
terns emerge when relatively simple systems are driven 
into unstable states that will deform dramatically in re- 
sponse to small perturbations. As the patterns are aris- 
ing from an instability, the pattern-forming behavior is 
likely to be extremely sensitive to small changes in sys- 
tem parameters. The description of deterministic pattern 
forming systems is typically accomplished in the form of 
partial differential equations such as the Navier-Stokes 
equations for fluids and reaction-diffusion equations for 
chemical systems. These phenomena have primarily been 
studied exclusively in continuum systems such as hydro- 
dynamical, optical, chemical systems and liquid crystals, 
although more recently pattern formation of similar type 
have also been reported in periodically vibrated granular 
media Q. 

Complementing the development of the theoretical un- 
derstanding of pattern formation in continuum systems, 
the localized mode forming ability of discrete lattices has 
also received significant recent attention. There is now a 
fairly complete understanding of the existence and sta- 
bility of localized structures (often referred to as intrinsic 
localized modes (ILM's) or discrete breathers) in a vari- 
ety of nonlinear lattices, undriven |^ as well as driven 
Q. It is fair to claim that these collective patterns are 
well understood while the process of their creation and 
interaction remains relatively unexplored. 

In the present communication we study the pattern 
forming abilities of a damped and periodically driven 
nonlinear lattice. Specifically, we demonstrate how a 
driven nonlinear (cubic) Klein-Gordon lattice forms a va- 
riety of patterns via modulational instabilities. We anal- 
yse the modulational instabilities and show how these re- 
late to the length-scale of the patterns that are formed. 
Normally, the spatial extend (characteristic length-scale) 
of ILM's is directly related to the frequency of the ILM's 
1^. However, in our case the (generally different) length- 
scale emerging from the instability may lead to a length- 
scale competion, the results of which we will explore. 



Model and stability of homogeneous solution. First, we 
study the Klein-Gordon lattice 



Xn + IXn + t^o^n = A„a;„ -f Aa;^ + e cos uot, 



(1) 



where 7 is the damping parameter, loq the natural fre- 
quency of the oscillators, A = ±1 the nonlinearity pa- 
rameter and finally e is the amplitude of the ac-drive at 
frequency oj. In one dimension the nearest neighbor cou- 
pling is AnXn = Xn+i ~ 2a;„ + Xn-i Q. The amplitude 
Aq (and phase i5o) of the spatially homogeneous solution 
Xn — Aocosiojt -\- 5q) of Eq.(|l|) can, within the rotating 
wave approximation, be shown to satisfy 



^0 7 



(2) 



The amplitude Aq of the response to a driving amplitude 
e is shown for the soft (A = 1) potential in Fig. 1. 




FIG. 1. Amplitude Ao of response vs driving amplitude 
e of the ac-drive at frequency, u = 1.2ajo (solid line) and 
uo = O.Sljo (dashed line). Remaining parameters are uo = 1-3, 
7 = 0.15a;o, and A = 1 (soft potential). 
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For u! < ujQ (dashed hne) three solutions are possible, 
while for oj > ujq (solid line) a single solution is possible. 
A similar picture is valid for the hard (A = —1) potential 
except that the multiple solutions then appear in the case 
w > Wo- 

Analyzing the stability of the homogeneous solution 
with respect to spatial perturbations, we introduce Xn = 
y + Zn into Eq.(|l|). Assuming periodic boundary con- 
ditions, we may expand z„ in its Fourier components 
Zn — exp(jfcn)^fc(i), where the mode amplitude ^fe(t) 
is then governed by 



£,k + i£.k + t^fcCfc 



-A^l^ [1 + cos(2wi + 25o)] ik (3) 



with — tjQ+4sin^(fc/2) denoting the linear dispersion 
relation of the system. 

Finally, the transformation ^fc(<) = C,k{'-^t + 
(5o) exp(— ^(wt + 5)) = Ckij) exp(--^r) reduces Eq.(||) to 
a standard Mathieu equation 



dr2 



+ aCfe-2gcos(2T)Cfe = 0, 



(4) 



where 



1 3A^2 
a= — {ALol~&\Al~^^) and g = (5) 



4^2 



As is well-known Q the Mathieu equations exhibit para- 
metric resonances when ^/a ~ z, where i = l,2,3,.... The 
width of the resonance regions depends on the ratio q/a 
(see, e.g. 0). In the framework of Eq.(^ the extent of 
the primary resonance a ~ 1 can easily M be estimated 
to be (a — 1)2 < . However, in the presence of the 
damping 7 the resonance condition for Eq.(|^) becomes 



q^>^,+{a-l)\ 



(6) 



With a and q defined in Eqs.(|^), given A, 7, 0^,0^0, and 
e, this translates into an instability band of certain 
wavenumbers k. 

Figure H shows this instability band as given by Eq. (|^) 
for parameters corresponding to the solid curve in Fig. 
|l|. The shaded region is the band of wavenumbers that 
are unstable according to Eq.(|6|) and the dashed line in- 
dicates the most unstable wavenumber, i.e. a = \. The 
effect of the damping clearly is to pinch off the instabil- 
ity region at a finite driving e > 0. Similar, using Eq.(^, 
instability regions can be determined for the solutions 
indicated by the dashed curve in Fig.|l]. 

We have verified the presence and location of the insta- 
bility band by direct simulations of Eq. (^ starting from 
the initial condition a;„(t = 0) = Aq cos{5q) + r/„, where, 
rjn, represents a small {\r)n\ ^ ^0) spatially random per- 
turbation. This initial condition injects energy into all 
wavenumbers and in the presence of an unstable region 
of wavenumbers the dynamics will enhance the energy 



content in this region and thereby identify the unstable 
region. 
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FIG. 2. Instability region (shaded) in the (e,fc)-plane. 
Dashed line indicates the most unstable wavenumber a = 1. 
Parameters corresponds to the solid line in Fig. 1. 

The above analysis is easily extended into the case 
of two spatial dimensions the only required change 
being that the dispersion relation now is Loi = lo"^ -\- 

4sin^(fcj,/2) -t- 4 sin^(A;y/2), where the wavevector is fc = 
(fcj;, ky). The instability in this case appears on an annu- 
lus in the wavevector plane, with a radius given by a = 1 
(see, Eq.(||)) and a width determined by Eq.(^). 

Pattern formation. Numerical simulations allow us not 
only to verify the predicted instability band, but also 
to follow the full nonlinear development and saturation 
initiated by the instability. In particular, in regions of 
parameter space we obtain the spontaneous formation 
of patterns of distinct spatial geometry. Although we 
have observed this phenomenon in one as well as in two 
dimensions, in the present communication we focus on 
the two-dimensional system, where the pattern formation 
is very rich. 

Although the dynamics show different features accord- 
ing to the specific region of parameter space, it is possible 
to trace a typical behavior as follows: Initializing the sys- 
tem in the spatially homogeneous state described above 
with a small amount of randomness added, the instability 
sets in after a certain number of cycles of the ac-drive, de- 
pending on the strength of the parametric resonance, i.e. 
on the value of ^/a ~ 1,2,3,.... Thereafter the system 
usually evolves through a sequence of different patterns 
(rhombi, stripes, etc.), composed of localized regions of 
high amplitude oscillations, before reaching a final con- 
figuration that may or may not result in a structure of 
definite symmetry. 

Due to the sensitive response to very small changes of 
the parameters, determining stability regions for the dif- 
ferent pattern geometries is a difficult task. However, in 
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the case of a hard potential (A = —1), Fig. 3 shows a hm- 
ited area of (e,a;)-space in which distinct spatial patterns 
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FIG. 3. Stability diagram showing the stability boundaries 
between the possible patterns in the (£,a;)-plane. Parameters 
are uo ~ 0.75, A = — 1 (hard potential), 7 = O.OSwo- The four 
points refer to the specific spatial structures shown in Fig. 4. 

This diagram is constructed following the full dynam- 
ics of the system for thousands of cycles. As our study 
concentrates on patterns arising from instabilities of the 
homogeneous solutions, we do not discuss possible hys- 
teretic behavior as is sometimes observed in similar sys- 
tems (see, e.g. Ref.[2]). 

Figure 4 shows representative examples of the spon- 
taneously emerging patterns corresponding to the points 
marked in Fig. 3. 



FIG. 4. (Color) Spatial patterns corresponding to the four 
points marked in Fig. 3 for e = 0.57: (a) straight stripes 
{u! = 0.8a;o), (b) modulated stripes {lo = O.SStJo), (c) rhombi 
{uj = LOwo), (d) localized rhombi {uj = 1.2ci;o). Xn.m is plotted 
along the vertical axes. 



The patterns consist of localized regions of high ampli- 
tude oscillations (ILM's) residing on a background that 
oscillates at the frequency, uj, of the ac-drive. In all the 
considered cases we have observed the natural result that 
patterns are only energetically sustained when the ILM, 
uJiLM, and driving, w, frequencies are commensurate, i.e. 
i^iLM = nuj, where n is an integer. For the patterns dis- 
played in Fig. 4, n — 2. Further, the motion of the 
ILM's is out of phase, i.e. , at the points in time where 
the background oscillation reaches its maximal excursion 
the ILM's obtain their minimal amplitude such that at 
these points the state is completely homogeneous. 

At a fixed driving e, for increasing values of the fre- 
quency to, as in Fig. 4, we can observe the following 
behavior: Due to the presence of the damping 7, at 
sufficiently small ui the spatially homogeneous solution 
is stable towards all possible spatial modulations such 
that the flat state is sustained. However, for values near 
point (a) the system becomes unstable with respect to 
certain spatial modulations and spatial patterns in the 
form of large stripes emerge (Fig. 4(a)). Increasing the 
frequency, these stripes become thinner and denser and 
begin to show an increasingly clearer modulation (Fig. 
4(b)). The characteristic length scale of these patterns is 
set by the size of the unstable fc- vector according to the 
above analysis. The nonlinear character of the system 
results in a transition towards a more isotropic geometry 
(rhombic) as the driving frequency is increased further 
(Fig. 4(c)). As in the case of the stripes, stronger lo- 
calization of the ILM's arranged in the rhombic pattern 
(Fig. 4(d)) is observed for even larger driving frequen- 
cies. The angle between the sides of the rhombus unit 
cell varies but for the hard potential it is always close to 
tt/2. For values of e larger than those displayed in Fig. 3, 
the final mesoscopic patterns of the system dynamics are 
spatially disordered much like the phenomena observed 
in granular media Q| . It is important to realize that the 
length of the unstable wavevector determines the length- 
scales of the final patterns, while the symmetry of the 
patterns is determined by the nonlinear character of the 
system. 

As a result of the periodic boundary conditions, the 
length scale of the emerging patterns must be commen- 
surate with the system size. However, we have observed 
that the existence of a band (variability in length and 
angle) of unstable fc-vectors (see Fig. 2 and related dis- 
cussion) allows continuous accommodation of this con- 
straint except for the discontinuous changes in length 
scales occurring when it is energetically favorable for the 
system to add (or subtract) an additional stripe (or row 
of ILM's). 

For the soft potential (A = 1) the variation of the am- 
plitude e, and the frequency, oj of the ac-drive is par- 
ticularly problematic as the dynamics in this case can 
lead to the development of catastrophic instabilities as 
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one or several oscillators overcome the finite barrier in 
the quartic potential.In all the cases we have been able 
to simulate, the early stage time evolution of the system 
is characterized by the formation of ILM's regularly ar- 
ranged in a square pattern. This spatial configuration, 
which is sustained for up to hundreds of cycles, seems 
always to suffer from a weak instability and eventually 
deforms into a rhombic pattern, as shown in Fig. 5. 



FIG. 5. (Color) Snapshot of spontaneously formed pattern 
of ILM's. The red line represents a one-dimensional cut of 
the pattern array, whose dynamics we have analyzed in detail 
(see Fig. 6). 

Contrasting with the case of the hard potential, here in 
the soft potential the angle between the sides of the rhom- 
bus unit cell is always close to 27r/3 (so almost hexago- 
nal). 

We now analyze this pattern more closely. The pat- 
tern shown in Fig. ^ consists of ILM's spontaneously 
organized into a regular rhombic pattern. The ILM's 
are spatially localized and perform harmonic temporal 
oscillations (at the frequency, w, of the ac-drive) and 
are therefore objects that are well described 1^ in the 
literature. A particular feature of these ILM's is that 
they reside on a background that oscillates at the same 
frequency. To expose the dynamics of the ILM's more 
closely, we show in Fig. ^ a time sequence along a one- 
dimensional cut of the two-dimensional system (the cut 
is indicated by the line in Fig. ||). It should be noted 
that in Fig. ^ we have removed the oscillations of the 
background in order to expose the ILM dynamics most 
clearly. 

Although we have focused on the ILM dynamics in the 
case of a soft potential, the features are analogous in the 
case of the hard potential and only the symmetry of the 
patterns is different. 



FIG. 6. ILM dynamics monitored along the direction indi- 
cated by the line in Fig. 5. 

In summary we have studied the modulational insta- 
bility in the damped and ac-driven cubic nonlinear Klein- 
Gordon lattice. The analysis applies to one as well as two 
spatial dimensions. We have further demonstrated how 
these instabilities lead to a variety of mesoscopic patterns 
of intrinsic local modes. In the case of a hard potential we 
characterized the patterns in a stability diagram and in 
the case of the soft potential we showed that the dynam- 
ics always result in a rhombic pattern with an angle close 
to 27r/3. These rhombic patterns were never observed in 
the case of the hard potential. This difference in the 
shapes of the rhombic patterns in the hard and soft cases 
can be understood by exploiting the analogy between the 
changes in the steady states of a dissipative systems and 
phase transitions in systems at thermodynamic equilib- 
rium (see. e.g. Ref. (H). In terms of this analogy the 
appearance of spatially periodic structures in a nonequi- 
librium system can be connected to a perturbation of 
the translational symmetry of thermodynamic states in 
equilibrium. A study based on this philosophy was pur- 
sued in Ref. with the results that the angle defining 
rhombic patterns is given by the coefRcient of the cubic 
term in the system, which is precisely our observation. 
From the present analysis it appears that experimental 
studies of the pattern forming abilities of discrete systems 
present an excellent opportunity to study ILM's and their 
mutual interactions and mesoscopic patterning. For ex- 
ample, optical systems |^ and spin systems po| appear 
good candidates for such studies. 

Research at Los Alamos National Laboratory is per- 
formed under the auspice of the US DoE. 
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